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Abstract 

Magnetostrictive  transducers  are  becoming  increas¬ 
ingly  prevalent  in  industrial  applications  including  high 
speed  milling  and  hybrid  motor  design  due  to  their 
broadband,  large  force  capabilities.  To  achieve  the  level 
of  performance  required  by  such  applications,  however, 
these  transducers  must  operate  in  nonlinear  and  hys- 
teretic  regimes.  To  accommodate  this  nonlinear  behav¬ 
ior,  models  and  control  laws  must  incorporate  known 
physics  and  be  sufficiently  robust  to  operate  under  re¬ 
alistic  operating  conditions.  We  develop  here  an  Hoc 
robust  control  design  for  a  prototypical  magnetostric¬ 
tive  transducer  but  note  that  the  design  is  sufficiently 
general  to  be  utilized  for  several  commonly  used  smart 
materials  including  piezoceramics  and  shape  memory 
alloys.  The  performance  of  the  control  strategy  is  illus¬ 
trated  through  a  numerical  example. 

1.  Introduction 

Applications  utilizing  piezoceramic  (PZT),  magne¬ 
tostrictive  and  shape  memory  alloy  (SMA)  compounds 
range  from  nanopositioning  stages  in  an  atomic  force 
microscope  (AFM)  to  vibration  suppression  systems  in 
buildings.  However,  the  material  attributes  which  pro¬ 
vide  these  compounds  with  unique  control  capabilities 
also  produce  hysteresis  and  constitutive  nonlinearities 
which  must  be  accommodated  in  models  and  control 
algorithms  to  meet  the  stringent  design  criteria  associ¬ 
ated  with  these  applications. 

The  method  we  employ  to  attenuate  this  nonlinear 
and  hysteretic  behavior  is  to  utilize  an  approximate 
inverse  of  the  hysteresis  model  as  a  filter  to  mitigate 
these  effects.  This  control  strategy  requires  a  hysteresis 
model  which  admits  an  inverse  that  can  be  employed  in 
real-time.  Filters  of  this  type  have  been  constructed  us¬ 
ing  a  domain  wall  model  [6,8]  as  well  as  Preisach  mod¬ 
els  [12].  Here  we  utilize  a  free-energy  based  hysteresis 


model  [9]  and  a  corresponding  approximate  inverse  al¬ 
gorithm  which  lends  itself  to  real-time  implementation. 

To  focus  the  discussion,  we  consider  an  industrial 
application  requiring  high  speed,  high  accuracy  milling 
using  the  prototypical  magnetostrictive  actuator  illus¬ 
trated  in  Figure  1.  The  objective  of  this  application 
is  milling  out-of-round  objects  at  speeds  of  3000  rpm 
within  a  tolerances  of  1-2  microns.  Stresses  and  dis¬ 
placements  are  provided  by  the  Terfenol-D  rod  in  re¬ 
sponse  to  applied  fields  generated  by  the  surrounding 
solenoid.  A  permanent  magnet  encasing  the  transducer 
provides  a  biased  field  to  achieve  bidirectional  strains 
as  well  as  providing  a  mechanism  for  flux  shaping. 

As  detailed  in  [1,  2],  these  transducers  possess  the 
ability  to  generate  broadband,  high  force  responses. 
However,  they  also  exhibit  hysteresis  and  nonlineari¬ 
ties  in  the  relation  between  the  input  field  H  and  the 
magnetization  M  generated  in  the  Terfenol-D  rod.  We 
focus  on  this  magnetostrictive  transducer  in  this  paper 
but  note  that  the  hysteresis  model  and  control  design 
are  sufficiently  general  to  permit  direct  extension  to 
analogous  piezoelectric  and  shape  memory  alloy  mod¬ 
els  as  developed  in  [3,  10,  11]. 

Whereas  employing  an  inverse  filter  attenuates  the 
primary  effects  of  the  nonlinear  and  hysteretic  behav¬ 
ior  of  the  smart  materials,  there  are  still  disturbances 
in  the  process  due  to  discretization  and  modeling  errors 
as  denoted  by  d  in  Figure  2.  In  realistic  applications, 
the  errors  d  are  not  the  only  disturbance  to  the  control 
system  which  must  be  accommodated.  Sensor  noise  is 
also  present  in  all  applications  and  can  cause  a  degra¬ 
dation  in  the  transducer  performance  if  unaccommo¬ 
dated.  Therefore,  we  consider  an  Hoo  design  capable 
of  rejecting  both  classes  of  disturbances. 


Figure  2.  Approximate  model  inverse  employed  as  a 
filter  for  robust  control  design. 


2.  Hysteresis  Model 

We  summarize  here  the  relevant  aspects  of  the  hys¬ 
teresis  model  presented  in  [9[.  This  model  quantifies 
the  energy  required  to  reorient  moments  in  combina- 
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tion  with  stochastic  homogenization  techniques  to  ac¬ 
commodate  variations  in  coercive  and  effective  fields. 
We  assume  fixed  temperatures  for  this  development  and 
the  model  ignores  losses  due  to  eddy  currents  and  there¬ 
fore  should  be  employed  for  operating  regimes  in  which 
eddy  current  losses  are  negligible.  A  further  assump¬ 
tion  is  that  the  spring  washer  in  the  transducer  provides 
sufficient  prestress  to  dominate  crystalline  anisotropies. 

We  first  quantify  the  internal  energy  due  to  the  inter¬ 
action  of  moments  through  the  Helmholtz  energy.  As 
detailed  in  [9] ,  under  the  assumption  that  moments  ori¬ 
ent  either  with  the  applied  field  or  diametrically  oppo¬ 
site  to  it,  a  reasonable  form  of  the  Helmholtz  energy  is 

r  i{m+mry 

^=J  l{M-MRy  ,  M  >  Mj  (!) 

I  ,  \M\  <  Mj. 

As  depicted  in  Figure  3,  Mr  and  Mi  respectively  de¬ 
note  the  point  at  which  the  minimum  of  ip  occurs  and 
the  inflection  point. 

In  the  presence  of  an  applied  field  H ,  the  Gibbs  en¬ 
ergy  is  given  by 

g  =  ip(M)  -  HM.  (2) 

We  note  the  magnetostatic  energy  is  £  =  Ho  HM,  where 
Ho  denotes  the  magnetic  permeability,  so  the  Gibbs  re¬ 
lation  (2)  can  be  interpreted  as  incorporating  no  into  A- 

Magnetoelastic  coupling  can  be  incorporated  through 
the  extended  Helmholtz  relation 

ipe{M,  e)  =  ip{M)  +  \yme2  -  YM'yeM2 

and  the  corresponding  Gibbs  energy 

G(H ,  M,  e)  =  ip(M)  +  i YM£ 2  -  Ym^eM2  -  HM  -  ae. 

Here,  ip  is  specified  by  (1),  Y M  is  the  Young’s  modulus 
at  constant  magnetization,  7  is  a  magnetoelastic  cou¬ 
pling  coefficient,  a  is  the  applied  stress  and  e  denotes 
the  strains. 

For  operating  regimes  in  which  thermally  activated 
processes  are  negligible,  the  local  average  magnetiza¬ 
tion  can  be  quantified  from  the  necessary  condition 
=  0.  In  the  absence  of  applied  stresses,  the  lo¬ 
cal  model  thus  predicts  a  linear  relationship  between 
H  and  M  with  slope  A  As  detailed  in  [9],  the  local 
magnetization  in  this  limiting  case  is  given  by 

r  0) 

[M(H,-Hc,0m  =  \  Y~Mr  (3) 

l  f  +mR 

for  the  respective  cases  {r(t)  =  0},  {t(£)  A  0  and 
i?(maxr(£))  =  —  Hc},  {r(£)  A  0  and  if  (max  r(t))  = 
Hc}.  The  transition  points  are  specified  as 

Tit)  =  {t  G  (0,T>]  I  Hit)  =  —Hc  or  Hit)  =  Hc) 


Figure  3.  (a)  Helmholtz  energy  ip  and  Gibbs  energy 
g  for  increasing  fields  H.  (b)  Dependence  of  the  local 
magnetization  M  on  the  field  H  at  the  lattice  level  in 
the  absence  of  thermal  activation. 


and  the  initial  moment  orientation  is 

r  f -MR,HiO)<-Hc 
[MiH-Hc,Om=  <  £  ,~Hc<HiO)<Hc 
{  f+MR,HiO)>Hc. 

For  a  description  of  models  which  incorporate  thermal 
activation,  the  reader  is  referred  to  [9]. 

The  relation  (3)  is  derived  under  the  assumption  that 
the  lattice  structure  is  homogeneous  which  is  overly 
simplistic  in  that  it  neglects  material  defects,  polycrys¬ 
tallinity  and  material  nonuniformities.  Also,  the  model 
assumes  that  the  effective  field  He  at  the  domain  level 
is  the  applied  field  H.  To  incorporate  these  effects, 
stochastic  distributions  are  employed  to  develop  a  bulk 
magnetization  model  for  nonhomogeneous  Terfenol-D 
samples  with  nonconstant  effective  fields  (see  [9]  for  de¬ 
tails). 

To  include  the  effects  of  material  nonhomogeneities, 
we  consider  a  distribution  of  free  energy  profiles  to  ac¬ 
commodate  the  nonhomogeneous  lattice  structures.  We 
assume  that  the  local  coercive  field  Hc  =  r]iMR  —  Mi) 
is  normally  distributed  with  mean  Hc.  Hence  the  co¬ 
ercive  field  has  the  density 

f(Hc)  = 

where  the  parameters  Ci  and  b  are  positive. 

The  second  extension  entails  the  incorporation  of 
variations  in  the  effective  field  through  the  assumption 
that  it  is  normally  distributed  about  the  applied  field 
H.  Therefore,  the  effective  field  has  the  density 

f(He)  =  C2e~^H-H ')2/5 

where  C2  and  b  are  positive.  Combination  of  the  coer¬ 
cive  and  effective  field  distributions  yields  the  magne- 
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tization  relation 


4.  Transducer  Model 


M(H)  =C  /  [ M(H ;  Hc,  £)]f(Hc)f(He)dHedHc. 
J  0  J  —  oo 

(4) 

The  integrals  are  approximated  using  a  compos¬ 
ite  quadrature  rule  for  the  numerical  implementation 
of  the  free  energy  model  thus  allowing  the  hysteresis 
model  to  be  computed  algebraically  and  hence  at  low 
cost.  Because  of  the  exponential  decay  of  the  distri¬ 
butions,  they  can  be  truncated  to  speed  computation. 
The  model  is  thus  sufficiently  efficient  to  allow  the  pos¬ 
sibility  of  real  time  implementation. 

The  elastic  constitutive  relation  is  provided  by  the 
equilibrium  condition  =  0  which  yields 

a  =  YM£-YM~/M2.  (5) 

The  coupled  constitutive  relations  (4)  and  (5)  quantify 
the  magnetization  and  stresses  for  an  undamped  mag- 
netostrictive  material.  These  relations  are  employed 
in  Section  4  when  quantifying  the  displacements  pro¬ 
duced  by  the  magnetostrictive  transducer  depicted  in 
Figure  1. 


3.  Inverse  Hysteresis  Model 

The  monotonicity  and  efficiency  of  the  hysteresis 
model  are  exploited  to  construct  an  approximate  in¬ 
verse.  To  determine  the  applied  field  required  to  create 
a  desired  magnetization,  the  hysteresis  model  is  used  to 
increment  the  magnetization  until  the  desired  magneti¬ 
zation  is  surpassed.  Then  the  applied  field  is  computed 
by  linear  interpolation.  The  computational  speed  of  the 
inverse  compensator  depends  on  the  size  of  the  step  AH 
taken  in  advancing  the  hysteresis  model.  Larger  steps 
will  increase  the  speed  while  decreasing  the  accuracy 
of  the  inverse  compensator.  The  relative  linearization 
error  for  an  input  signal  with  a  frequency  of  1  Hz  and 
a  step  size  of  AH  =  1  employed  in  the  inverse  model  is 
plotted  in  Figure  4.  A  successful  control  design  must 
be  able  to  reject  this  error  to  the  input  of  the  plant. 
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Figure  4.  Relative  error  created  by  inverse  filtering. 


In  this  section,  we  employ  the  relations  (4)  and  (5)  to 
quantify  the  displacements  generated  by  applied  fields 
in  the  prototypical  transducer  depicted  in  Figure  1. 
The  model  for  the  transducer  is  based  on  the  model 
developed  in  [2]. 

One  end  of  the  rod  (x  =  0)  is  assumed  to  be  fixed 
while  the  other  end  ( x  =  L )  is  constrained  by  a  damped 
oscillator  and  has  an  attached  point  mass  as  depicted  in 
Figure  5.  The  internal  damping  coefficient  and  density 
of  the  Terfenol-D  rod  are  denoted  by  cM  and  p  and  the 
point  mass  is  denoted  M l.  The  end  spring  has  stiffness 
Ul  and  damping  coefficient  Cl- 

The  constitutive  relation  (5)  is  for  an  undamped 
material.  To  incorporate  Kevin-Voigt  damping,  the 
stresses  at  any  point  x,  0  <  x  <  L,  are  assumed  to  be 
proportional  to  a  linear  combination  of  strain,  strain 
rate  and  squared  magnetization  which  yields 

a  =  YM^+cMp^_  yM  M  2 
ox  oxot 

where  the  linear  strain  relation  £  =  §y  is  employed. 

Balancing  forces  yields 


82w  =  3Ntot 
dt 2  dx 


(6) 


where  A  is  the  cross  sectional  area  of  the  Terfenol-D 
rod  and  the  force  resultant  Ntot  =  JA  cdA  is  given  by 


Ntot 


YmA 


dw 

dw 


cM  A 


d2w 

dxdt 


YmA^M2. 


To  obtain  appropriate  boundary  conditions,  we  first 
note  that  w(t,  0)  =  0.  As  detailed  in  [2],  balancing 
forces  at  x  =  L  gives 


n  rj2 

Ntot(t,  L)  =  -kLw{t,L)  -  cL-^(t,L)  -  MLg^-(t,L). 


Initial  conditions  are  taken  to  be  w(0,x)  =  0  and 

^(0,*)=0. 

The  model  can  be  implemented  utilizing  the  weak 
form  and  a  Galerkin  finite  element  approximation.  To 
further  simplify  our  system,  we  assert  that  the  magne¬ 
tization  in  the  Terfenol-D  rod  can  be  taken  as  uniform 


+W 


x=L 


^/VvV^ 

— H — 

p  ifw 
n.  d  t 


Figure  5.  Boundary  conditions  for  the  rod  model. 
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Figure  6.  System  representation  including  input  disturbances  d  and  sensor  noise  n. 


over  the  length  of  the  rod.  This  is  reasonable  in  present 
actuator  designs  since  flux  shaping  via  the  surround¬ 
ing  magnet  is  used  to  minimize  end  effects  in  the  rod. 
Hence,  since  each  section  of  the  rod  reacts  identically  to 
the  uniform  magnetic  field,  the  elements  act  uniformly 
and  therefore  the  transducer  dynamics  can  be  mod¬ 
eled  as  a  damped  spring-mass.  To  achieve  bidirectional 
strains,  we  linearize  the  magnetostrictive  relationship 
about  a  biasing  magnetization  level  of  Ms/2.  This  bias 
can  be  achieved  my  means  of  the  permanent  magnet  as 
depicted  in  Figure  1.  For  more  details  on  this  model 
development  see  [2,  5]. 

The  dynamics  of  the  Terfenol-D  transducer  can  be 
now  represented  by 

x  +  kx  +  cx  =  uM(t) 

(7) 

a;(0)  =  Xq  ,  i(0)  =  Xq. 


The  signal  n  represents  noise  in  the  sensor  the  mea¬ 
surements  of  y.  For  the  simulation  results  presented  in 
Section  6,  we  assume  sensor  noise  with  a  frequency  of 
60  Hz.  The  output  signals  e  and  u  denote  the  weighted 
tracking  error  and  weighted  output  of  the  controller 
K,  respectively.  The  weighting  functions  Wu,  Wd ,  We, 
Wr  and  Wn  are  chosen  to  maximize  the  performance 
of  the  controller  utilizing  a  priori  knowledge  regarding 
the  characteristics  of  the  signals. 

We  now  construct  the  transfer  function  representa¬ 
tion  of  the  open-loop  system.  The  transfer  matrix  from 
the  inputs  r,  d ,  n  and  u  to  the  outputs  e,  u  and  v  is 
specified  by 


G  = 


WeWr  -WePWd  -WeWn  -WeP 

0  0  0  wu 

Wr  -PWd  —Wn  -P 


The  scalars  k,  c,  and  v  are  determined  by  fitting  the 
model  (7)  to  the  Galerkin  approximation  of  (6)  or  to 
data  from  the  physical  device.  For  our  sample,  k  = 
7.8899  x  103,  c  =  6.4251  x  107  and  v  =  1.3724  x  1CT2 
yielded  accurate  model  fits.  Note  that  the  hysteresis 
inherent  to  the  Terfenol-D  rod  is  still  present  in  (7) 
through  (4)  which  quantifies  the  hysteretic  relationship 
between  H  and  M. 


The  control  system  can  then  be  represented  as  a  linear 
fractional  transformation  as  shown  in  Figure  7.  Details 
regarding  the  system  formulation  can  be  found  in  [5,  7]. 
To  proceed,  we  partition  G  as 


A 

Bi 

b2 

G(a)  = 

Ci 

0 

P>12 

c2 

P>2\ 

0 

5.  Robust  Control  Design  where 


In  this  section,  we  present  a  robust  control  design  for 
the  Terfenol-D  transducer  depicted  in  Figure  1.  The 
control  system  incorporates  the  presence  of  external 
disturbances,  such  as  the  errors  cause  by  the  inverse 
filter  and  sensor  noise,  and  minimizes  their  effects  with 
respect  to  the  Hoo  norm 


G  n  = 

G21  = 


A 

Bi 

Ci 

0 

A 

Bi 

c2 

D2i 

A 

b2 

Ci 

P>ii 

A 

b2 

c2 

0 

(9) 


Moo  =  sup  a[T(jw)] 

ue  r 

where  a[T(ju>)\  denotes  the  maximum  singular  values 
of  the  closed-loop  map  T.  Analogous  H2  control  laws 
are  developed  in  [5,  7]. 

Figure  6  illustrates  the  block  diagram  of  the  system 
to  be  controlled.  In  the  diagram,  P  represents  the 
transducer  model  given  by  the  differential  equation  (7). 
The  signal  to  be  tracked  and  the  position  of  the  tip  of 
the  Terfenol-D  rod  are  respectively  denoted  by  r  and 
y.  The  signal  d  represents  the  error  in  the  lineariza¬ 
tion  of  the  input  by  the  inverse  filter  (see  Figure  4). 


respectively  represent  the  transfer  functions  from  w  to 
z,  u  to  z,  w  to  v,  and  u  to  v,  (see  Figure  7).  To  en¬ 
sure  the  existence  of  a  sub-optimal  controller,  the 
following  assumptions  are  made  regarding  the  system: 


1.  (A,  Hi)  is  controllable  and  (Ci,  A)  is  observable, 

2.  (A,  B2)  is  stabilizable  and  (C2,  A)  is  detectable, 

3.  D*12Di2  >  0  and  D2iD^  >  0, 

A  —  jujl  B2 


4. 


Ci 


D, 


has  full  column  rank  for  all  w, 
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z= 


=  w 


Figure  7.  Linear  fractional  transformation  (LFT)  rep¬ 
resentation  of  transducer. 


5. 


A  —  jujl  B\ 

C-2  D21 


has  full  row  rank  for  all  to. 


Details  regarding  the  validity  of  these  assumptions  are 
provided  in  [5,  7]. 


5.1.  Weighting  Functions 

The  selection  of  the  weighting  functions  is  critical 
to  the  performance  of  the  robust  control  design.  A 
discussion  of  the  choice  of  these  function  can  be  found 
in  [4,  5,  7,  13].  By  design,  the  frequency  of  the  reference 
signal  r  is  taken  to  be  1  Hz  so  we  construct  the  pass- 
band  filter  Wr  to  have  a  bandwidth  of  1  Hz  centered  at 
1  Hz.  The  frequency  of  the  noise  n  can  be  accurately 
determined  for  the  device  measuring  the  position  of  the 
tip  of  the  Terfenol-D  rod.  For  the  simulation  results, 
noise  with  a  frequency  of  60  Hz  is  added  to  the  system 
so  Wn  is  taken  to  be  a  pass-band  filter  with  a  bandwidth 
of  10  Hz  centered  at  60  Hz. 

To  determine  the  weighting  function  Wd,  a  signal 
with  the  same  frequency  as  the  reference  signal  was 
filtered  by  the  inverse  compensator  and  then  fed  into 
the  hysteresis  model.  The  output  is  defined  to  be  the 
desired  control  signal  plus  the  disturbance  d  (see  Fig¬ 
ure  2).  The  power  spectrum  of  d  indicates  that  the 
most  significant  frequencies  in  d  lie  below  400  Hz  so 
Wd  is  taken  as  a  low-pass  filter  with  a  cut-off  frequency 
of  400  Hz. 

The  weighting  function  on  the  error  signal  was  taken 
to  be  We  =  with  qe  =  4  x  106  and  ee  =  1  x  10-8. 
An  integrator  was  chosen  to  prevent  the  error  from 
achieving  steady  state  at  a  nonzero  value  and  the  pole 
was  shifted  slightly  off  zero  to  ensure  that  the  controller 
design  is  realizable.  We  specified  the  weight  on  the  con¬ 
troller  output  to  be  Wu  =  5  x  10-6.  Since  we  do  not 
experience  any  problems  with  saturation  or  other  such 
effects,  we  minimally  weight  u  to  focus  the  controller  on 
tracking  and  disturbance  rejection.  It  is  important  to 
note  that  increasing  the  order  of  the  filters  increases  the 
number  of  states  in  the  controller  design.  To  facilitate 
real-time  implementation,  a  minimal  control  realization 
was  utilized  to  limit  the  number  of  states  in  the  control 
gain  K. 


5.2.  Sub-Optimal  Control  Design 

Employing  the  notation  defined  in  (9),  the  design  of 
a  sub-optimal  controller  which  gives  ||T||oo  <  7 
incorporates  two  Riccati  equations 

A* X  +  XA  +  X  -  B2B*2)  X  +  C*C i  =  0 

(10) 

and 

AY+YA*+Y  (w2C*Cn  -  C*C2 )  Y+B1Bl  =  0.  (11) 

The  following  theorem  from  [13]  guarantees  the  exis¬ 
tence  of  an  TLoo  sub-optimal  control. 

Theorem  1:  There  exists  an  admissible  controller  such 
that  llTH-o  <  7  if  and  only  if 

1.  There  exists  a  positive  definite  solution  X  to  (10), 

2.  There  exists  a  positive  definite  solution  Y  to  (11), 

3.  p{XY)  <  72. 

The  TLoo  optimal  controller  is  subsequently  given  by 


A 

-ZL 

F 

0 

with 


A  =  A  +  7  ~2BxBIX  +  B2F  +  ZLC2, 
F=~B*2X,  L=-YC |. 

The  control  given  in  (12)  is  sub-optimal  in  that  it 
provides  a  closed-loop  system  with  an  Hoo-norm  less 
than  7.  The  hinfsyn  command  in  Matlab  can  be  uti¬ 
lized  to  decrease  7  until  an  assumption  of  Theorem  1 
is  violated.  This  yields  a  control  design  method  which 
is  close  to  optimal. 

5.3.  Numerical  Example 

We  present  here  a  numerical  example  illustrating  the 
Hoo  robust  control  law.  The  results  were  computed 
using  a  noise  signal  n  with  a  magnitude  of  1  x  10— 5 , 
which  is  10%  of  the  reference  signal,  and  a  frequency  of 
60  Hz.  The  inverse  compensator  was  computed  using  a 
step  size  of  AIL  =  1  in  the  approximate  inverse  model 
to  ensure  a  simulation  of  the  control  process  which  has 
the  potential  for  real-time  implementation.  Figure  8a 
illustrates  the  controller’s  ability  to  track  the  reference 
signal  and  reject  the  noise  and  disturbance  signals.  A 
tracking  error  less  that  2  microns  is  achieved  after  a 
short  period,  as  illustrated  in  Figure  8b.  Additional 
examples  illustrating  the  performance  of  the  method 
are  provided  in  [7]. 
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Figure  8.  (a)  TLoo  tracking  performance  and  (b)  tracking  error. 


6.  Concluding  Remarks 

An  Hoq  robust  control  design  for  a  magnetostrictive 
system  has  been  summarized  in  this  paper.  It  is  shown 
that  the  control  design  is  capable  of  maintaining 
accurate  tracking  while  rejecting  sensor  noise  and  a  dis¬ 
turbance  due  to  an  inexact  inverse  filter.  The  same 
methodology  can  be  used  to  include  other  disturbances 
to  the  system  if  necessary.  Finally,  while  the  methods 
were  illustrated  in  the  context  of  a  magnetostrictive 
transducer,  they  are  sufficiently  general  to  be  applied  to 
systems  utilizing  piezoceramic  or  shape  memory  alloys 
due  to  the  unified  modeling  framework  used  to  quan¬ 
tify  hysteresis  and  constitutive  nonlinearities  inherent 
to  all  of  these  materials  [10]. 
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